Radial variation of attenuation and star formation in the largest 
late-type disks observed with GALEX 
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ABSTRACT 

For a sample of 43 nearby, late-type galaxies, we have investigated the radial variation of both 
the current star formation rate and the dust-induced UV light attenuation. To do this we have 
cross-correlated IRAS images and GALEX observations for each of these galaxies, and compiled 
observations of the gas (CO and HI) and metal-abundance gradients found in the literature. 
We find that attenuation correlates with metallicity. We then use the UV profiles, corrected 
for attenuation, to study several variants of the Schmidt law and conclude that our results are 
compatible with a simple law similar to the one of Kennicutt (1998a), extending smoothly to 
lower surface densities, but with considerable scatter. We do not detect an abrupt break in the 
UV light at the threshold radius derived from Ha data (at which the Ha profile shows a break 
and beyond which only a few HII regions are usually found). We interpret the Ha sudden break 
not as a change in the star formation regime (as often suggested) but as the vanishingly small 
number of ionizing stars corresponding to low levels of star formation. 

Subject headings: galaxies: spiral, ultraviolet: galaxies, infrared: galaxies, (ISM:) dust, extinction 

1. Introduction 

Interstellar dust affects our view of galaxies 
from the UV to the near-infrared. High-redshift 
galaxies are commonly studied in the rest-frame 
UV where dust effects can be extremely severe, 
and our estimation of the "cosmic" star formation 
rate (SFR) (e.g. Schiminovich et al. 2005) crucially 
depends on the corrections applied to account for 
the dust attenuation of starlight. 

Radiative transfer models suggest that the ratio 
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of far-infrared to UV radiation is a reliable mea- 
sure of the UV attenuation (Buat & Xu 1996; Witt 
& Gordon 2000; Panuzzo et al. 2003), depending 
weakly upon the geometry of stars and dust, the 
extinction law, or the nature of the underlying 
stellar population. Since this ratio is not always 
available, the slope of the UV spectrum redden- 
ing with respect to local starbursts (Calzetti et 
al. 1994; Meurer et al. 1995, 1999) has been com- 
monly used as a metric to estimate the amount 
of attenuation. This allows one to compare the 
attenuation (as measured by the far infrared/UV 
ratio) to the slope of the UV spectrum in a di- 
agnostic plot, the commonly-called "IRX-/3 rela- 
tionship" . Several recent works, mostly based on 
GALEX observations, cast doubt on the IRX-/3 re- 
lation that has been most commonly used during 
the recent years (Buat et al. 2005; Cortese et al. 
2006; Scibcrt et al. 2005; Gil dc Paz et al. 2006). 
Independent studies have come to similar conclu- 
sions, e.g. Bell (2002) for normal spirals and Bell 
et al. (2002) for individual regions in the LMC. 

These works suggests either that there is no re- 
lationship (or a very noisy one) or that a rela- 
tionship exists but differs from the starbursts one, 
depending on the selection of the sample (infrared 
vs UV or optically-selected). 

The largest angular sized galaxies in the 
GALEX Atlas of Nearby Galaxies (Gil de Paz 
et al. 2006) allow us to revisit the IRX-/3 rela- 
tionship. Indeed, several tens of nearby galaxies 
observed by GALEX are large enough to have been 
resolved even by IRAS. The combination of the 
UV data of GALEX and infrared of IRAS opens 
up the possibility of looking at the spatially re- 
solved IRX-/3 relationship across a relatively large 
sample of galaxies. Even if the number of resolved 
objects in this study is smaller than the total 
number of objects used in the global studies, each 
of our galaxies provides a number of independent 
points along the radial profile. Regions situated 
at the same radius within a galaxy are relatively 
homogeneous (for instance similar metallicities), 
thus radial profiles are physically relevant. 

Another compelling reason to study attenua- 
tion^ in these galaxies is that, as our closest neigh- 

"Attenuation" is the loss of flux due to the presence of dust. 

It is sometimes considered that "extinction" should be used 
only for the dust properties and not the global effect of dust 



bors, they have been extensively surveyed at other 
wavelengths, and in the published literature there 
are metallicity gradient and gas density profiles, 
etc. that allow us to go one step further and try 
to characterize how the properties of dust might 
depend on these other physical quantities. We also 
produce star formation rate radial profiles derived 
from the UV (and not Ha data), corrected for the 
dust attenuation for all our galaxies. These pro- 
files can be tiscd to study how the star formation 
rate depends on the gas density, i. e. the so-called 
Schmidt law. We will be able to use these star for- 
mation rate profiles to finally study the threshold 
radius beyond which little (if any) star formation 
is observed. Many studies suggest that the thresh- 
old radius corresponds to the point where gas den- 
sity becomes lower than a critical density for star 
formation (e.g. Martin & Kennicutt 2001). It is 
important to bring as many constraints as possi- 
ble to bear on the complex problem of star forma- 
tion and this paper brings its share, by using UV 
GALEX data. Indeed, while many other studies of 
the Schmidt law (Kennicutt 1998a; Boissier et al. 
2003; Wong & Blitz 2002) are based on Ha data, 
the star formation rate in this paper is derived 
from the UV, which "averages" the star forma- 
tion on longer time-scales than Ha. Recombina- 
tion lines as Ha are emitted only in presence of 
a large ionizing flux, coming mostly from massive 
stars (M > IOMq) with lifetimes shorter than 2 
X lO*" years. On the contrary, the UV continuum 
is emitted by stars with lifetimes up to 10* years 
(a review of the various star formation indicators 
and their associated lifetimes is given in Kennicutt 
1998b). Because of the scarcity of UV data before 
GALEX, little work had been done in the UV on 
the star formation law. Some exceptions arc, that 
of Buat et al. (1989), Bersier et al. (1994) con- 
cerning respectively M81 and M51, or Buat (1989) 
and Donas et al. (1987) for studies using sample 
of nearby galaxies but these papers looked only 
at the integrated properties of galaxies, and em- 
ployed rather simple recipes for dust correction. 
The radial variation of the UV-derived star for- 
mation rate, correctly dust-corrected via the far- 
infrared to UV ratio, is the subject of this paper. 
This is an important work as we now know that 



on the galaxy light, also depending on e.g. the geometry 
of the system. "Extinction" is nevertheless often used with 
the same meaning as "attenuation" . 
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the UV also shows us low levels of star formation 
(XUV galaxies, e.g. Gil de Paz et al. 2005; Thilker 
et al. 2005), in regions with such low dust con- 
tent that the far-infrared emission is undetected 
by IRAS, and that the UV is likely to be weakly 
affected by extinction. 

In section 2. we discuss how the data-sets were 
assembled and prepared in order to undertake this 
study of the radial variation of the attenuation, of 
the star formation, and of their relation to other 
physical quantities in those same galaxies. In sec- 
tion 3, we present our study of the attenuation de- 
rived from these data, and the consequences con- 
cerning the resulting star formation law are pre- 
sented in section 4. 
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Table 1 

Sample: basic parameters used and references 



Name 


E(B-V) 


T 


D 


PA 




(mag) 




(Mpc) 


(deg) 


IC 1613 


0.025 


10 


0.8 


50 


IC 2574 


0.036 


9 


4.0 


50 


M31 


0.062 


3 


0.8 


35 


M33 


0.042 


6 


0.8 


23 


M51a 


0.035 


4 


8.4 


-17 


M64 


0.041 


2 


17.0 


-65 


M81 


0.080 


2 


3.6 


-23 


M83 


0.066 


5 


4.5 





MlOl 


0.009 


6 


7.5 


90 


NGC 0055 


0.013 


9 


2.0 


-72 


NGC 0247 


0.018 


7 


3.1 


-6 


NGC 0253 


0.019 


5 


3.9 


52 


NGC 0300 


0.013 


7 


2.0 


-69 


NGC 0628 


0.070 


5 


11.2 


25 


NGC 0660 


0.065 


1 


12.4 


-10 


NGC 0891 


0.065 


3 


9.6 


22 


NGC 0925 


0.076 


7 


9.3 


-78 


NGC 1097 


0.027 


3 


15.3 


-50 


NGC 1291 


0.013 





9.7 


-15 


NGC 1512 


0.011 


1 


10.4 


90 


NGC 1566 


0.009 


4 


17.5 


60 


NGC 2366 


0.036 


10 


3.4 


25 


NGC 2403 


0.040 


6 


3.2 


-53 


NGC 2841 


0.016 


3 


14.1 


-33 


NGC 2903 


0.031 


4 


8.9 


17 


NGC 3198 


0.012 


5 


16.7 


35 


NGC 3521 


0.058 


4 


9.0 


-17 


NGC 3621 


0.080 


7 


8.3 


-21 


NGC 3627 


0.032 


3 


9.1 


-07 


NGC 4216 


0.032 


3 


17.0 


19 


NGC 4236 


0.015 


8 


4.5 


-18 


NGC 4395 


0.017 


9 


4.2 


-33 


NGC 4438 


0.028 





17.0 


27 


NGC 4559 


0.018 


6 


17.0 


-30 


NGC 4569 


0.046 


2 


17.0 


23 


NGC 4594 


0.051 


1 


9.1 


90 


NGC 4631 


0.017 


7 


9.0 


86 


NGC 4656 


0.013 


9 


9.0 


33 


NGC 4736 


0.018 


2 


5.2 


-75 


NGC 5055 


0.018 


4 


8.2 


-75 


NGC 7331 


0.091 


7 


14.9 


-09 


NGC 7793 


0.019 


7 


2.0 


-82 



D25 d25 Res. References 



(') 


(') 


(") 


[0/H] 


HI 


CO 


V(R) 


16 


15 


106 


d 


24 




24 


13 


5 


110 


e 


25 




25 


190 


60 


108 


i 


28 


28 


18, 27 


71 


42 


108 


i 


17 


19 


17 


11 


7 


117 


i 


5 


30 


40 


10 


5 


116 




10 


36 


10 


27 


14 


111 


i 


5 


5 


40 


13 


12 


106 


i 


6 


6 


6 


29 


27 


88 


b 


23 


23 


23 


32 


6 


96 


i 


33 




33 


21 


7 


107 




14 




14 


28 


7 


114 


i 


34 


41 


34 


22 


16 


106 


i 


32 




32 


11 


10 


111 


a, h, i 


39 


36 


39 


8 


3 


108 






46 


43 


14 


3 


111 




38 


37 


37, 38 


11 


6 


110 


h, i 


5 


5 


31 


9 


6 


103 


g 




46 


29 


10 


8 


88 




42 






9 


6 


95 










8 


7 


95 


i 




2 




8 


3 


95 


f 


20 




20 


22 


12 


103 


h, i 


5 


5 


3 


08 


4 


105 




3,9 


45 


3 


13 


6 


106 


h, i 


5 


5 


3 


9 


3 


103 


i 


4 




8 


11 


5 


104 


i 




36 


15 


12 


7 


105 


i 








9 


4 


110 




47 


36 


47 


8 


2 


108 




16 


22 




22 


7 


102 




11 






13 


11 


104 


h,f 


44 




44 


9 


3 


106 










11 


4 


104 


i 


12 




12 


10 


4 


110 




7 


7 


7 


9 


4 


108 




1 




1 


16 


3 


106 








35 


15 


3 


106 










11 


9 


104 


i 


26 


36 


26 


13 


7 


105 


i 


9 


36 


9 


11 


4 


104 


i 


3 


45 


3 


9 


6 


108 


i 


13 




13 
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Table 1 — Continued 
WLM 0.037 10 1.0 -20 12 4 104 c 21 



Note. — The columns give (1) Object name, (2) Galactic color excess (Schlegel et al. 1998) , (3) RC3 
morphological type, (4) Distance (Mpc) taken from the GALEX Atlas of nearby galaxies (Gil de Paz et al. 
2006), (5) Position Angle (degree), (6) D25: major axis diameter (minutes), (7) d25: minor axis diameter 
(minutes), (8) FWHM of the restored 100 microns beam in arcsec (see section 2 for details), (9) References 
for 0/H gradients, (10), (11), (12) References for HI, CO, and rotation curve data. 

References. — References for Abundance Gradients [O/H]: a BcUcy & Roy (1992) ; b Kcnnicutt 
et al. (2003a) ; c Lee et al. (2003) ; d Lee et al. (2003) ; e Masegosa et al. (1991) ; f Roy et al. (1996) ; g 
Storchi-Bergmann et al. (1996) ; h van Zee et al. (1998) ; i Zaritsky et al. (1994). References for Gas [HI, 
CO] and Kinematic [V(R)] Data: 1 Bajaja ct al. (1984) ; 2 Bajaja ct al. (1995) ; 3 Bcgcman (1987) ; 
4 Begeman (1989) ; 5 Boissier et al. (2003), and references within ; 6 Boissier et al. (2005), and references 
within ; 7 Boselli et al. (2006), and references within ; 8 Bosma (1978) ; 9 Bosma (1981) ; 10 Braun et al. 
(1994) ; 11 Braun (1997) ; 12 Brocils & van Woerden (1994) ; 13 Carignan & Puchc (1990a); 14 Carignan 
& Puche (1990b) ; 15 Casertano & van Gorkom (1991) ; 16 Cayatte et al. (1994) ; 17 Corbelh (2003) ; 18 
Dame et al. (1993) ; 19 Heyer et al. (2004) ; 20 Hunter et al. (2001) ; 21 Jackson et al. (2004) ; 22 Kenney 
& Young (1988) ; 23 Kenney et al. (1991) ; 24 Lake & Skillman (1989) ; 25 Martimbcau ct al. (1994) ; 26 
Mulder & van Driel (1993) 27 Newton & Emerson (1977) ; 28 Nieten et al. (2006) ; 39 Ondrechen et al. 
(1989) ; 30 Paglione et al. (2001) ; 31 Pisano et al. (1998) ; 32 Puche et al. (1990) ; 33 Puche, Carignan & 
Wainscoat (1991) ; 34 Puchc, Carignan & van Gorkom (1991) ; 35 Rand (1994) ; 36 Regan ct al. (2001) ; 
37 Sakamoto et al. (1997) ; 38 Sancisi & Allen (1979) ; 39 Shostak & van der Kruit (1984) ; 40 Sofue (1997) 
; 43 Sorai et al. (2000) ; 44 van Driel et al. (1988) ; 45 van Driel et al. (1995) ; 46 Wevers et al. (1986) ; 48 
Young & Scoville (1982) ; 50 Young et al. (1995) ; 51 Zhang et al. (1993). 
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2. Profiles: Data and Procedure 

Our initial sample consists of 48 late-type 
galaxies (RC3 T type going from to +10, thus 
including all spirals from SO/a to Sm and magel- 
lanic irregulars Im) with UV data in the GALEX 
Atlas of nearby galaxies (Gil do Paz et al. 2006), 
larger than 8 arcmin (as measured by D25 diam- 
eter corresponding to the B-band isophote of 25 
mag arcsec"^.) This size allows for spatial res- 
olution of these galaxies at all the wavelengths 
needed, including the 60 and 100 /^m IRAS data. 

All the IRAS images were obtained using 
HIRES requests to the IPAC web page^. The 
method used in HIRES is described in Aumann, 
Fowler and Melnyk (1990); examples of images ob- 
tained with HIRES have been presented by Rice 
(1993) for a sample of nearby galaxies (The HIRES 
procedure failed to converge for NGC 3109, and 
this galaxy was removed from our analysis). Be- 
cause of the IRAS aperture itself and the complex 
pattern of observations, the resulting HIRES im- 
ages can be quite asymmetric. HIRES produces 
simulated images of point sources distributed into 
the image at the sky position of the targeted ob- 
ject. The first step in our procedure was to apply 
an asymmetric Gaussian convolution in order to 
produce a circular image at 60 /im and 100 /im of 
these point sources, and to apply the same con- 
volution to the HIRES reconstructed image of the 
galaxy. 

In the UV, the GALEX FUV and NUV (Far 
and Near UV, respectively around 1516 and 2267 
A) images of the Atlas (Gil de Paz ct al. 2006) were 
used. For four galaxies, no FUV data were avail- 
able (NGC 1365, NGC 3628, M98, M106), and 
they are excluded from our analysis which in fine 
includes 43 galaxies. Given this selection process, 
this sample is not complete in any sense, although 
it should be representative of optically selected 
galaxies (a comparison of the GALEX Atlas of 
nearby galaxies and the Nearby Field Galaxy Sur- 
vey of Jansen et al. 2000 is given in Gil de Paz et 
al. 2006). 

Stars were removed using the star masks of 
the Atlas and then interpolated over. Foreground 
Galactic extinction was corrected as in Gil de Paz 
et al. (2006) : using the values of the color excess 

^http: / /irsa.ipac. caltech.edu/IRASdocs/hires_over. html 



E(B-V) given by the maps of Schlegel et al. (1998) 
and a Galactic extinction curve. 

Images in FUV, NUV, and at 60 /xm were then 
convolved with an elliptical Gaussian function in 
order to match the resolution to the one mea- 
sured in the HIRES simulation of point-sources 
at 100 /im (after the asymmetric Gaussian convo- 
lution discussed above). The spatial resolution of 
IRAS is relatively poor. However, this might help 
us to avoid diSiculties in interpreting the UV/IR 
balance since the dust of small regions could be 
heated by the UV emission from neighboring re- 
gions. 

Profiles were finally computed within elliptical 
annuli with ELLIPSE in IRAF*^ using convolved 
images, with fixed center, ellipticity and position 
angles (as given in Gil de Paz et al. 2006 and 
Table 1). The step in major axis radius between 
isophotes was chosen to be the measured resolu- 
tion (close to 1.5 arcmin on average, and given in 
Table 1). 

The infrared profiles were combined as in Dale 

et al. (2001) in order to compute a "far-infrared" 
profile {FpiR in W m~^ arcsec"^) and a "total- 
infrared" profile (Ftir), as described below: 

Ffir = 1.26 X 10-14(2.58/60 + /loo) (1) 

{fx are the IRAS flux surface densities in Jy 
arcsec"^) and 

logiFrm) =log{FFiR) -F0.2738 - 0.0282r 0.7281r2 

+0.6208r^ + 0.9118r* (2) 

where r = log{fm/ fim) 

Attenuation profiles (expressed in magnitude) 
were then computed using the fit of Buat et al. 
(2005) 

A{FUV) = -0.0333X3 _^ 0.3522^2 

+1.1960X + 0.4967 (3) 

where X = log{FTiR/ Fpuv) ; Fpuv = ^fuvfi^f^^ 
and fv;^^ is the FUV flux surface density in W 
m~^ Hz~^ arcsec"^. Other expressions linking 
A{FUV) and X do exist (e.g. Burgarella et al. 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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2005) with subtle differences, and since this kind 
of caUbration might also be subject to future revi- 
sions, we included in our figures both scales (ob- 
served Ftir/FfuVi calibrated A{FUV)). 

Given the long chain of operations in comput- 
ing A{FUV), we estimated in a simple way its un- 
certainty by computing the extremal values it can 
take when moving the UV and infrared fluxes by 
1 a from their observed values. FIR and UV data 
are moreover uncertain by ~ 15 % on average, (see 
e.g. Morrissey et al. 2005 for an assessment of the 
quality of the GALEX photometry, and Rice 1993 
for the photometry of IRAS extended galaxies). 

The additional uncertainty in the calibration of 
the relation given in equation 3 itself is not taken 
into account. Points and error-bars were not al- 
lowed to include negative attenuations. In the few 
cases where the A{FUV) obtained in this way was 
unphysical (e.g. for numerical reason), the last 
consistent estimated error was used for points at 
larger radii (this only affects a few outer points 
where error-bars are anyway large as can be seen 
in the A{FUV) profiles given in the appendix). 

The FUV observed profiles were finally cor- 
rected for attenuation using A{FUV), corrected 
for inclination, and converted to a star formation 
rate using the Kennicutt (1998b) calibration. The 
inclination correction accounts for the fact that 
the thickness of a galaxy (and thus the amount 
of stars) along the line of sight has to be multi- 
plied by cos{inc) (inc being the inclination of the 
galaxy) to obtain the same values for a face-on 
orientation (mc=0). The adopted inclination is 
given by the ratio of D25 and d25 from Table 1, as- 
suming the objects arc intrinsically circular. Note 
that this procedure assumes that we computed the 
correct A{FUV) for the inclination of the galaxy. 
Panuzzo ct al. (2003) showed that the relation be- 
tween A{FUV) and the Fpm/ Fpuv ratio actu- 
ally depends on inclination. Their fit for edge-on 
galaxies differs from the fit for moderately inclined 
galaxies by up to about 0.3 magnitude of extinc- 
tion at 2000 A. There are other reasons why we 
should be cautious about very inclined galaxies: 
i) Our data points are separated by one element 
of resolution along the major axis. In the case 
of inclined galaxies, the ellipses along the minor 
axis will be closer to each other than the resolu- 
tion. Outer regions might therefore be contami- 
nated by inner regions, ii) In an edge-on disk, the 



amount of material along the line of sight is such 
that only a full radiative-transfer modeling can be 
used to reconstruct the UV to sub-millimeter spec- 
tral energy distribution of galaxies (e.g. Popescu 
et al. 2000). For this reason, our data are pre- 
sented when possible with different symbols for 
various axis ratio ranges, and we discussed in the 
text when the inclination affected our results. 

Corollary data were collected from the litera- 
ture, concerning the gaseous content (HI, CO, Ro- 
tation curve) and the oxygen abundance gradient 
of the galaxies. CO data were converted to molec- 
ular H2 using the conversion factor of Boselli et 
al. (2002), dependent on the metallicity when an 
abundance gradient was available. For the seven 
galaxies without metallicity gradient, we used the 
conversion factor corresponding to the H-band lu- 
minosity of each galaxy, using the calibration of 
Boselli et al. (2002). 

The gas profiles taken from the literature were 
converted to the same inclination as mentioned be- 
fore (indicated by the index "inc"). We also deter- 
mined face on values (index "0" ) of the same quan- 
tities performing a similar correction than the one 
applied to the UV surface brightness (equivalent to 
a star formation rate surface density). Note that 
face-on values should be used to compare star for- 
mation rate and gas surface densities; while the ex- 
tinction profiles obtained for the observed galaxy 
inclination should be compared to the column den- 
sities along the line of sight, i.e. the ones with the 
"inc" index. While we corrected the original gas 
data to have consistent inclination with our own 
UV and far- infrared computations, we must ac- 
knowledge that this compilation introduces some 
level of uncertainty since all observations were not 
performed with the same instruments, spatial res- 
olution, etc. We however consider that this ef- 
fect is included within the typical uncertainties 
discussed below. 

Typical imcertainties are 0.15 dex for abun- 
dances, about 10 % for HI data and 50 % for H2. 
The fact that we used a constant CO-to-H2 con- 
version factor when we had no abundance gradi- 
ent for 7 galaxies must introduce an error in their 
molecular gas column densities. However, Fig. 1 
of Boissier et al. (2003) shows that usually the 
difference between a constant and a metallicity- 
dependent conversion factor is lower than 0.5 dex. 
When computing the total gas density (HI-I-H2), 
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a significant number of points have only HI data, 
and the contribution of H2 was then assumed to 
be negligible. This of course introduces another 
uncertainty on the total gas density, however: i) 
Most of those points correspond to outer regions of 
disks in which CO is detected and measured only 
for the inner portion of that same galaxy. It is 
reasonable to assume that the molecular fraction 
is small in these outer regions, ii) In 13 galaxies, 
we only have HI data (even in the center of the 
galaxy) but taking into account that only a minor 
fraction of the total gas is usually in the form of 
molecules (about 15 %, sec Boselli ct al. 2002) and 
that when we do have both profile, H2 dominates 
only for a few central points, and in most cases 
by only a factor a few; we must conclude that the 
error we commit is relatively small in most cases. 
We will nevertheless keep in mind these drawbacks 
when discussing our results. 

When we found several sets of data, we usu- 
ally adopted one (that seemed better -more re- 
cent, better sensitivity), and checked the consis- 
tency between them. In a few cases, we adopted 
an average or a combination of several references. 
Note that we interpolated the profiles given in the 
literature at the resolution adopted in our work. 
The original resolutions are very variable, but usu- 
ally better than the IRAS one. The gaseous pro- 
files obtained following this procedure are given 
in the figures in the appendix for each individual 
galaxy (at the same resolution as for the rest of 
our work). Table 1 gives the references for all of 
these data. 

3. Attenuation in the disk of late type 
galaxies 

3.1. Attenuation radial profiles 

The attenuation A{FUV) profiles are shown for 
each individual galaxy in the appendix. Most of 
them present a global decrease of extinction with 
radius. At large radii, error-bars are very large 
due to low levels of measured infra-red emission, 
so that the data are always consistent with no at- 
tenuation in the outer disks. The central value 
of the estimated FUV attenuation is typically of a 
few magnitudes (ranging from virtually to about 
6 magnitudes). The most central point of our pro- 
files could be affected by the presence of a bulge 
with e.g. UV upturn, or of a AGN. The fact that 



we sec in the individual A(FUV) profiles continu- 
ous trends with radius let us think that on average 
our extinction gradients are not much affected by 
such cases (in the figure of the appendix, only M31 
and NGC0253 do present a strong disturbance in 
the central part of the A{FUV) profile). 

A similar trend of decreasing extinction gra- 
dients was present in the six galaxies of the 
FOCA/IRAS study of Boissier et al. (2004). Hol- 
werda et al. (2005b) used counts of distant galax- 
ies, with a "synthetic field method" (Holwerda et 
al. 2005a) to estimate extinction gradients. Their 
totally independent results are consistent with 
ours, although the uncertainties are larger in their 
method. 

At this point, we would like to include some 
warning concerning the derived extinction A{FUV) 
By adopting the Buat et al. (2005) calibration, 
we assume that the star formation histories they 
adopted are representative of real galaxies. How- 
ever, it is probably not the case for a few galax- 
ies (or at least, parts of galaxies) suffering heavy 
quenching of the star formation rate (that could be 
induced, for instance by ram pressure in clusters, 
e.g. Boselli et al. 2006). In such galaxies, the 
GALEX FUV-NUV color would be intrinsically 
red (close to or above 1), with a different relation- 
ship between Ftir/ Fpuv and A{FUV) than the 
one we adopted. Such situations will be studied in 
Cortese et al. (in preparation). Part of the scatter 
in the right part of Fig. 1 could be due to such 
effects. As a result, applying extinction correction 
to individual galaxies (or profiles) should still be 
done with caution, especially for red FUV-NUV 
colors. In our approach, the number of galaxies 
is sufficiently large so that such effects should be 
minor on our global results. 

3.2. The Infrared Excess - UV Slope 
(IRX-/3) Relationship 

Heckman et al. (1995) and Meurer et al. (1995, 
1999) demonstrated with lUE data the existence 
of a relation between the slope of the UV spec- 
trum {13) and the TIR / UV ratio (so called IRX- 
(3 relationship) for starburst galaxies. It has been 
common thereafter (especially for high redshift ob- 
jects) to estimate the extinction from the UV slope 
using the IRX-/3 relationship (e.g. Schiminovich 
et al. 2005). In the recent years, however, several 
studies have shown that this relation does not hold 
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for all galaxies. For example, Bell (2002) showed 
that it fails in nearby galaxies, and Buat et al. 

(2005) came to the same conclusion for samples 
of UV and far-infrarcd-sclcctcd galaxies. Kong et 
al. (2004) proposed models with a large variety of 
star formation histories, demonstrating that they 
could span a large portion of the IRX-/3 plane. 
In this framework, the starburst relation would 
hold for galaxies with extremely large current-to- 
past-average star-formation rates (birthrate b pa- 
rameter), while galaxies with less activity would 
lie systematically below the relation. However, it 
seems that a relation can be found in nearby galax- 
ies (see e.g. Gil de Paz et al. 2006; Cortese et al. 
2006; Seibert et al. 2005), given a simple offset 
with respect to the original relation proposed for 
starbursts. This means that the range of star for- 
mation histories among real galaxies is probably 
not as wide as among models of Kong et al. (2004) 
and that the b parameter might not be the driver 
of the IRX-/3 relationship. Actually, Burgarella et 
al. (2005) showed that the slope of the extinction 
and the presence or not of a bump has a larger 
impact on the position in the IRX-/3 plane than 
the star formation history (and the b parameter). 
As discussed in Gil de Paz et al. (2006), Cortese et 
al. (2006) and Seibert et al. (2005), the difference 
in the IRX-/3 relation for starbursts and normal 
galaxies might be due to geometrical effects, aper- 
ture effects present in the lUE data. The relative 
calibrations in the computation of /3 with lUE's 
data and GALEX FUV-NUV color index could 
also play a role. More work remains to be done. 

In this context, our nearby galaxies have the ad- 
vantage of being close enough so that we can spa- 
tially resolve the IRX-/3 relationship. While this 
has been done by looking at individual regions in a 
few galaxies already, for instance M51 (Calzetti et 
al. 2005), we present for the first time the IRX-/3 
relationship for the profiles of all the large late- 
type galaxies of the GALEX Nearby Galaxies At- 
las in Figure 1. It is particularly interesting to 
note that while outliers are present (with corre- 
spondingly larger error-bars), most of the points 
describe a very tight relation, but shifted from 
the classical relation derived for/from starbursts, 
similar to the recent works of Gil de Paz et al. 

(2006) ; Cortese et al. (2006); Seibert et al. (2005) 
(see discussion above for the possible reasons for 
this shift). In order to quantify the relation- 



ship, we performed a nonlinear least squares fit 
of the form y = log{W+^'' - c) (with a;=FUV- 
NUV, the color between the two galex bands, and 
y—log{FTiR,/ Fpijv), after rejecting quite edge- 
on galaxies(NGC 0055, NGC 0253, NGC 0660, 
NGC 0891, NGC 4216, NGC 4621, NGC 4656a). 
We obtained for a, b, c respectively 0.570, 0.671 
and 3.220 using y-uncertainties as weight, and 
0.561, 0.713, 3.136 using no weights. These fits 
and their associated 1 sigma dispersion are shown 
in the right part of Fig. 1. 

The fact that our IRX-/3 relationship is quite 
tight is also interesting as there is much more scat- 
ter when individual regions are looked at, see for 
instance Calzetti et al. (2005) in M51a, or Thilker 
et al. (2006) in NGC 7331, both studies based 
on the G ALEX/SINGS joint analysis of GALEX 
and Spitzer data from the SINGS project (Kenni- 
cutt et al. 2003b). By azimuthally averaging pro- 
files, we remove the effects due to small-scale star- 
formation histories (affecting the UV slope), and 
is likely to reduce the effect of radiative transfer 
peculiarities on small scales: that is, the dust in 
some small regions may well be heated by the UV 
emission from neighboring regions. Using the low 
resolution of IRAS (coupled with azimuthally av- 
eraging) may turn out to be an advantage in this 
respect. 

Figure 1 also presents the integrated values for 
the galaxies of the Atlas (Gil de Paz et al. 2006) as 
grey squares. The squares overploted with a cir- 
cle correspond to the integrated value given in the 
Atlas for the galaxies of the present study. The re- 
lation of our study seems tighter than the one for 
the integrated galaxies of the Atlas. This can be 
explained by the combination of several effects: i) 
For all our galaxies, we followed exactly the same 
procedure from scratch, i.e. we do not use pub- 
lished fluxes but rather re-derive them ab initio, 
in contrast the UV Atlas for which fluxes were 
compiled from a wide variety of sources, ii) In in- 
tegrated galaxies, centers with activity might af- 
fect the whole flux if they are strong enough, while 
they probably influence only the central point in 
a profile like the one we compile. About half of 
our galaxies have a sign of central activity (AGN, 
LINER, Seyfert indicated in NED), but precisely 
because it is confined to the central point, it should 
not affect much our profiles. 
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3.3. The Absence of a trend between At- 
tenuation and Gas Column Density 

It is relatively common to consider the atten- 
uation to gas density ratio as being constant. 
For instance, Komugi ct al. (2005) adopt Nh = 
2 X lO^^Ay. This is essentially inspired by stud- 
ies in the Milky Way such as Bohlin et al. (1978). 
Under this assumption, one would expect a strong 
correlation between gas density and attenuation, 
what is not supported by our results (see Fig. 2, 
3). Thus, attenuation (and therefore the existence 
of dust) is probably affected by various processes 
not directly related to the gas density. In Fig. 
2, we show the amount of extinction as a func- 
tion of the gas in the form of neutral hydrogen 
(HI, left) and in the form of molecular hydrogen 
(H2, right). Within individual galaxies (points 
connected by one dotted curve), a trend of in- 
creasing attenuation with gas densities is gener- 
ally visible (specially with H2). However, the dif- 
ferences from one galaxy to another are large, so 
that globally no relation is clearly emerging. The 
fact that the amount of attenuation docs not seem 
to correlate with the HI colum density was also 
found by Holwerda et al. (2005b). A small trend 
can be seen between the molecular gas and the 
extinction. While this seems reasonable (inner re- 
gions are denser, have larger extinction and higher 
molecular fraction) and is observed in most galax- 
ies individually, the variations from one galaxy to 
another are still extremely large. In Fig. 3, we 
compare the ultraviolet attenuation A{FUV) to 
the total gas density. In this case, we are in an 
intermediate situation. At low gas column densi- 
ties, which correspond usually to the outer part 
of galaxies (where errors are larger), no trend is 
observed. In the inner parts of galaxies, the situa- 
tion is the same as for H2 since the molecular gas 
often dominates the total gas density or at least 
appears to compensate for the central plateau or 
the observed decrease of HI in the central regions 
of galaxies. 

In conclusion, our data do not support any re- 
lation of the type A/Nh = constant (commonly 
used in the literature). This is actually not sur- 
prising as A{FUV) is not a measure of the dust 
mass, but of the amount of dust heated by nearby 
young massive stars. As such, it will be largely 
affected by geometrical effects, and the physical 
properties of the dust grains (type, size distribu- 



tion...), what would probably erase any underlying 
relation between the dust and gas masses. Because 
the dust mass is dominated by cold dust, A{FUV) 
is just not a good measure of it. Estimating the 
dust mass, and thus the dust to gas mass ratio is 
out of the scope of this paper since IRAS wave- 
lengths do not allow us to probe the cold dust. 

3.4. The Attenuation-Metallicity Rela- 
tionship 

While the relation observed with the gas is ill- 
defined, the attenuation seems to depend in a 
clearer way on the metallicity (Fig. 4). Indeed, 
for a given abundance of oxygen, the scatter in ob- 
served extinction is smaller than at any given gas 
surface density. This is especially true if we ex- 
clude the galaxies with the largest inclinations in 
which the method used to derive the attenuation 
might suffer from various problems as commented 
in section 2. We performed a least squares mini- 
mization fit to the points within galaxies with axis 
ratio larger than 0.4 and found: 

A{FUV) = 1.02 X (12 + log{0/H)) - 7.84. (4) 

(the dispersion around this fit is (7a{fuv)=^-^ 
mag, and the correlation coefficient is 0.65). 
This is qualitatively consistent with the earlier 
work of Boissier et al. (2004). Albeit a bit shal- 
lower, this relation is characterized by similar at- 
tenuations in the high metallicity range (at 12 -|- 
log{0/H) ~ 9.3). Due to the limited number of 
galaxies, there were no points at the lower metal- 
licities (7.5-8.5) in Boissier et al. (2004). A source 
of potential difference between this earlier study 
and the present one is that it was performed at 
the FOCA wavelength rather than GALEX. In 
Boissier ct al. (2004), the extinction determined in 
the integrated star forming galaxies of Buat et al. 
(2002) with metalHcities of Gavazzi et al. (2004) 
were also shown to exhibit a similar trend. Our 
results are also in good agreement with the recent 
study of Cortese et al. (2006) (although their trend 
is slightly stronger than ours), for a sample of in- 
tegrated cluster galaxies observed with GALEX 
(their fit of the TIR-to-FUV ratio is converted to 
our A{FUV) scale in Fig. 4). 

Note that the attemiation-metallicity relation- 
ship we found concerns galaxies forming stars in 
a relatively quiescent (disk)-mode. Heckman et 
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al. (1998) have shown that the UV slope is corre- 
lated with the metalUcity in starbursts, also cor- 
responding to a metallicity-extinction relationship 
among starburst, shifted however with respect to 
the spirals to higher extinctions (and presented in 
Boissier et al. 2004). 

4. Star Formation Law 

4.1. The Theoretical Star Formation Laws 

The star formation on galactic scales is obvi- 
ously a crucial phenomenon in the evolution of 
galaxies. As a result, it is also a fundamental el- 
ement of any galaxy evolutionary model, e.g. N- 
body /semi-analytical models like GALICS (Hat- 
ton et al. 2003), SPH/chemical evolution models 
(Lia et al. 2002), chemical / spectrophotometric 
evolution models (Boissier & Prantzos 1999). De- 
spite its importance, only rough theories exist and 
a few empirical relationships have been looked for 
(e.g. Madore 1977; Buat et al. 1989; Kennicutt 
1998a; Wong & Blitz 2002). An excellent review 
of the empirical situation is given in Kennicutt 
(1998b), while Elmegreen (2002) consider the pos- 
sible physical origin of the empirical laws. The 
reader is refereed to these works for an extensive 
discussion on the topic. Here, we only come back 
rather quickly on a few of the relations proposed 
for the "star formation law" . These "laws" should 
allow us to predict the star formation rate from 
other physical quantities, as shown below. Espe- 
cially, we will look for relations between the star 
formation rate surface density (Ssfk) and the 
gas surface density in its various phases (neutral, 
molecular) as well as in total gas (neutral -|- molec- 
ular), under the form of a traditional "Schmidt 
law" (see Schmidt, 1959 for the original work): 

HsFR = aS^AS- (5) 

It is reasonable to expect a more direct rela- 
tionship with the molecular gas rather than the 
neutral gas since this phase is more closely in 
the sequence related to the ensuing star formation 
event. For instance, in their disk model consisting 
of self-gravitating clouds interacting to produce 
a turbulent viscosity, Vollmer & Beckert (2002) 
found that the star formation rate and molecu- 
lar gas surface densities should follow the same 
dependence with galactic radius. On the obser- 
vational side, while a relation between star for- 



mation rate and molecular mass was not initially 
confirmed (Kennicutt 1989; Boselli et al. 1995), 
Boselli et al. (2002) showed that using a luminosity 
or mctallicity-dcpendcnt CO-to-H2 conversion fac- 
tor provides a better relation than using HI alone. 

While this suggests a more direct connection 
between the star formation rate and the molecu- 
lar gas, the cause for star formation could still be 
related to the total gas. Indeed, a basic model 
of self gravitating disk would suggest a Schmidt 
law with the total gas density as the main control 
parameter (Kennicutt 1998a). 

The other model we will consider here is a 
Schmidt law modulated by a dynamical factor: 

i:sFR = aJ:-aAS^ (6) 

where V{R) is the rotational velocity at radius R 
(where the surface densities are determined). The 
ratio V{R)/R can bo considered cither as repre- 
senting the frequency of the passage of density 
waves enhancing star formation (Ohnishi 1975; 
Wyse & Silk 1989) or a dynamical time-scale (Ken- 
nicutt 1998b). In their model of the chemical 
evolution of the Milky Way, Prantzos & Aubert 
(1995) found that an another factor in addition 
to a simple Schmidt law was necessary and opted 
for a relation of the type given by equation 6 with 
V{R) = a constant, and n — \. Subsequent mod- 
els of the Milky Way and spirals used a similar 
SFR law (Boissier & Prantzos 1999, 2000) but 
with n = 1.5. This was found to give agreement 
with the Ha profiles of 16 nearby spirals studied 
by Boissier et al. (2003) although other laws could 
not be excluded. Boissier et al. (2003), had also 
tested a dependence on the stellar surface density 
in addition to the gas, as proposed by Dopita & 
Ryder (1994). Despite having an additional free 
parameter, this formulation did not offer a bet- 
ter fit to their data. To explore this further we 
would need near-infrared surface brightnesses for 
our galaxies. Unfortunately 2MASS data are not 
sufficiently deep in outer parts of our galaxies, and 
we decided to omit consideration of this law in our 
present study. 

Before testing the above proposed laws, we 
should note that the field of galactic star formation 
is not lacking new ideas. For instance, Elmegreen 
(2005) proposed that a combination of processes 
trigger star formation, conforming to the observed 
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relations and thresholds (densities below which no 
star formation would occur). Seigar (2005) pre- 
sented an interesting correlation between the spe- 
cific star formation rate and the shear rate, stress- 
ing the possible role of shear in star formation; 
however the shear versus average star formation 
rate (his Figure 4) has large dispersion. Our rota- 
tion curves do not have the necessary resolution or 
homogeneity to explore trends involving the shear 
rate. It has also been recently suggested (Barnes 
2004) that star formation may depend not only on 
the local gas density but may also be induced by 
shocks, especially in the case of interacting galax- 
ies (we do not have such violently interacting sys- 
tems in our sample). Testing this kind of star for- 
mation law requires detailed dynamical modelling 
on a case by case basis. Barnes (2004) showed 
that this idea is promising by reproducing obser- 
vations of the Mice. It would be however hard to 
implement such a law in e.g. semi-analytic model 
since it needs a detailed dynamical modelling of 
each interacting systems (including their precise 
geometry and details of the encounter). 

In any case, a best possible law for star forma- 
tion should still properly reproduce the relation 
observed between gas, dynamics and star forma- 
tion indicators that we compiled for our GALEX 
galaxies, whatever the physical causes are. Thus, 
this purely empirical work provides additional con- 
straints for future theoretical work and galaxy evo- 
lution models. 

4.2. Our constraints on the star formation 
laws 

In each panel of Fig. 5 we show a variant of 
the simple Schmidt law (SFR versus neutral gas, 
molecular gas, total gas), testing equation 5 and a 
combination of the SFR + dynamical factor ver- 
sus total gas testing the relation of equation 6. In 
this figure, the internal contributions of individual 
galaxies are connected by dotted lines. Individ- 
ual SFR and gas profiles, as well as SFR vs gas 
are given for each galaxy individually in the Ap- 
pendix. Correlation coefficients are indicated in 
each panel, as well as the two (dashed) regression 
lines whose coefficients are given in Table 2. For 
each relation, the first regression line is obtained 
by performing a least square fit of the star for- 
mation rate surface densities with the gas surface 
densities being given (minizing the errors on the 



Y axis), and the other one by performing a least 
square fit of the gas surface densities with the star 
formation rate surface densities being given (min- 
imizing the errors on the X axis). 
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Table 2 

Parameters a, and bi (i=l,2) op the two regression lines {log{Y) = Ui log{X) + bi) shown as 

DASHED LINES IN EACH PANELS OF FIGURE 5 (LEAST SQUARE FIT MINIMIZING THE ERRORS ON X AND Y 

respectively). We also give for each set of variables {X,Y) the correlation coefficient r 

AND THE NUMBER OF POINTS USED FOR THE FIT {N). 



X 


Y 


ai 


bi 


a2 


62 


r 


N 






0.83 


-0.30 


3.04 


-1.11 


0.52 


232 




^SFR 


0.57 


0.58 


1.26 


0.63 


0.67 


108 




^SFR 


0.99 


-0.44 


2.09 


-0.92 


0.68 


236 




SsFflxR/ V(R) 


0.68 


-1.53 


1.76 


-2.04 


0.62 


228 
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It is obvious from this figure that the form of the 
observed star formation laws differs from one spi- 
ral galaxy to another. This was already known and 
discussed earlier (e.g., Boissier et al. 2003; Kenni- 
cutt 1998a). It is hoped that averaging over many 
galaxies is equivalent to averaging the radial star 
formation history of any given galaxy over time. 
The trend with H2 seems relatively good on the 
high density side (larger than ~ 1 M0 pc~^) but 
the number of points with molecular gas alone is 
very small, and the results are inconclusive, espe- 
cially given the additional uncertainty on the con- 
version factor. As in the case of extinction vs II2, 
however the scatter seems too large to be caused 
by this factor alone. We note that Komugi et al. 
(2005) show that II2 is very well correlated with 
the star formation rate but they are interested in 
high-density regions and galactic nuclei, working 
at high resolution. In our case, to the contrary, 
the central parts of the galaxies have little influ- 
ence in our plots since we are working at low res- 
olution and looking at the largest possible radii. 
When the total gas is used, as in the study of the 
attenuation, we acknowledge that the absence of 
molecular gas data can affect the precise value of 
the slope and the intercept; however, the scatter 
seems too large to be entirely due to the uncer- 
tainties in the II2 content alone. Boissier ct al. 
(2003) showed that the extrapolation (or not) of 
CO data, and the use of a metallicity-dependent 
or a constant conversion factor did change signifi- 
cantly their results, but mostly influenced the star 
formation threshold (without much affecting the 
rate exponent n). 

Among all of the Schmidt law variants that we 
tried, the simple Schmidt law (with a dependence 

on the total gas density) gives the best results, 
although the correlation coefficient for the molec- 
ular gas and the "dynamical" law (respectively 
0.67 and 0.62) arc close to the one for the total 
gas (0.68). Komugi et al. (2005) have shown that 
the simple Schmidt law extends to higher gas sur- 
face densities when smaller and central regions (as- 
sumed to be dominated by H2) are considered. If 
added to our diagrams, their data would go from 
about 10 to 1000 MqPC^'^ for the gas, correspond- 
ing to a range in star formation rate from about 
1 to 1000 Mgpc-^Gyr-^ This is compatible with 
an extrapolation of our adopted Schmidt law (us- 
ing the total gas) to larger surface densities, with 



a similar slope. We notice however a shift to lower 
star formation rates for the same gas density. In 
their study, the SFR is derived from Ha, which 
might explain this offset (see below for a discus- 
sion on this point). 

In all panels of Fig. 5, we also show for compar- 
ison the results of Boissier et al. (2003), in the form 
of a hourglass-shaped shaded area (the extent of 
the hourglass being the extent of the observations, 
and the diagonals the regression lines found in this 
work). This study used Ha profiles for 16 spiral 
galaxies. There are some systematic differences: 
on average they found lower star formation rates 
for the same gas density than in our present study, 
except with HI. The slope of a simple Schmidt law 
was also slightly steeper than what we find. Over- 
all, taking into account the large dispersion, the 
differences are not dramatic and may be caused by 
i) inadequate/uncertain corrections: for instance, 
in Boissier et al. (2003) Ha was corrected for ex- 
tinction and [Nil] contamination by standard fac- 
tors depending only on the galaxies type (a com- 
mon procedure), but not considering possible ra- 
dial variations, ii) calibration issues : the conver- 
sion of Ha and UV light to star formation rates 
depends on the initial mass function adopted for 
converting the observed fluxes (Kennicutt 1998b). 
For example, if the assumed IMF overestimates 
the number of ionizing stars, the star formation 
rate derived from Ha will be underestimated. An- 
other difference, that cannot be explained by these 
effects is the fact that the gas extent was smaller, 
i.e. with our UV profile, we detect low levels of 
star formation, corresponding to lower levels of gas 
surface densities. 

In Fig. 6, wc compare the simple Schmidt law 
(obtained from our UV radial profiles versus to- 
tal gas) to the galaxies of Kennicutt (1998a), ob- 
tained from Ha imaging of spiral galaxies, aver- 
aged within their optical disk. Considering the 
scatter, the agreement between both studies is 
quite good. Here again, however, we have many 
more points at low gas densities. We also in- 
clude in this diagram the star formation rate val- 
ues derived from the UV assuming zero attenua- 
tion (crosses). These are lower limits on the star 
formation rate, showing that a significant amount 
of star formation is detected at low gas densities, 
and this result does not depend on the attenuation 
correction. This comes from several facts: i) Ken- 
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nicutt's estimates arc averages within R25, where 
densities are on average higher ii) In the UV, we 
do find emission at low levels of star formation (in 
stark contrast to Ha observations), corresponding 
to gas surface densities below the putative thresh- 
old for star formation. 

In this figure, wc also connected the points cor- 
responding to M31. The inner 50" of this galaxy 
are orthogonal to the usual Schmidt relation (and 
the one found in the rest of the sample) . We leave 
for later studies to interpret this untypical behav- 
ior of one of our nearest neighbor. 

Our plots might appear more highly dispersed 
than other studies (e.g. Komugi et al. 2005; Kenni- 
cutt 1998a) but this is mainly because we do not 
work on a large range of gas densities (i.e., the 
scatter only seems smaller in Kennicutt (1998a) 
because of the larger dynamic range introduced 
into the plots when circumnuclear starbursts are 
included), and we are working primarily at low 
densities where stochastic effects and larger uncer- 
tainties may have a larger impact. Another pos- 
sible cause of the dispersion is that we are using 
UV data, and as Buat et al. (1989) noted, the 
time-scales we are probing is significantly longer 
than with Ha, probably making the underlying 
"causal" relations harder to catch (Madorc 1977). 

4.3. Concerning the Notion of a Threshold 

Star formation rate profiles, derived from Ha 
data are known to frequently present an abrupt 
break, while atomic gas is generally much more 
spread than the star forming and the stellar disc 
(e.g. Martin & Kennicutt 2001). This observa- 
tion is traditionally interpreted as the result of lo- 
cal instabilities developing in differentially rotat- 
ing discs within a threshold radius where the gas 
density is larger than a critical density (Toomre 
1964). More complex physics has also been sug- 
gested (e.g. the onset of thermal instability of 
Schaye 2004). 

The above-described observation of a Schmidt- 
like behavior at very low gas densities as derived 
from UV radial profiles suggests that low-levels 
of star formation are in fact common beyond the 
usual threshold radius, where Ha emission ceases 
or is not (as easily) observed. We can verify that 
this is actually the case by taking the 9 galaxies for 
which we have UV and infrared data. Ha profiles 



from Martin & Kennicutt (2001), and for which 
the threshold radius given in Martin & Kennicutt 
(2001) is larger than 90 arcsec. This last condition 
ensures that we would resolve the threshold at our 
low resolution. The star formation rate profiles de- 
rived from the UV and the Ha emissions are pre- 
sented in Figure 7. Within the threshold radius 
Rthreshoid, they are in rough agreement, while a 
systematic difference occurs beyond it: eight out 
of the nine galaxies show no truncation (or any 
change of slope) in the UV profiles at the posi- 
tion of the Ha threshold. The last one is M81, 
for which our UV profiles ends at the position of 
the threshold. Note that we are limited for this 
galaxy by its large angular size: the south outer 
part of the galaxy is getting close to the limit of the 
GALEX field of view. The infrared background 
is also relatively bad so that we cannot extend 
the profiles very far out in the disk. However, an 
inspection of the FUV image (at full resolution) 
clearly reveals two spiral arms extending in the 
area between R/ RthreshoU equals 1 and 1.5. 

The apparent absence of threshold may seem at 
odds with the figures presented in the appendix, 
in which an abrupt drop of the UV surface bright- 
ness (or SFR surface density) is visible for a few 
galaxies. This is the case for instance of NGC0891. 
However, for this galaxy, the UV drop is concomi- 
tant with a drop in the gas density, and occurs 
at quite large radius (35 kpc). Various reasons 
may affect the gas at large radii (e.g. gravita- 
tional interactions, ram-pressure), and obviously 
no star formation would occur when the gas reser- 
voir was removed. Thus, while abrupt drops ex- 
ist in the UV for some galaxies, they do not al- 
ways correspond to a "threshold" of star forma- 
tion. NGC0925 presents a more interesting case 
where no UV is detected beyond about 20 kpc, 
while HI is still observed. Even if this indicates 
for this galaxy a "UV threshold" , our analysis 
includes galaxies with detected star formation at 
similar gas densities. Detailed studies of individ- 
ual galaxies with high rc;solution may help to un- 
derstand the peculiarities of system like this. This 
kind of analysis is beyond the scope of this paper. 

Despite these questions on possible "UV thresh- 
olds", it is clear from figure 7 that beyond the 
usual threshold radius, we frequently see UV radi- 
ation when Ha ceased to be steadily found. While 
many reasons for this are possible, including a 
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variation in the initial mass function, the simplest 
explanation is that at some radius, star forma- 
tion levels are so low that the number of high- 
mass (short-lived) ionizing stars visible at any time 
becomes vanishingly small, while slightly longer- 
lived stars, emitting in the UV for a correspond- 
ingly longer time, arc more frequent. If this in- 
terpretation is correct, the origin of the thresh- 
old of star formation, largely debated on theoret- 
ical grounds (Toomre 1964; Quirk 1972; Wang & 
Silk 1994; Martin & Kennicutt 2001; Wong & Blitz 
2002; Boissier et al. 2003; Schaye 2004) on the ba- 
sis of Ha observations would be simply the lack of 
ionizing stars. 

To test further this idea, we took the Ha profiles 
of Martin & Kennicutt (2001) around the thresh- 
old radius they defined. Since profiles are usually 
computed within annular ellipses, we converted 
them to star formation rate integrated whithin 1 
kpc wide annular ellipses (that we will call tpi), 
what is shown in figure 8. For galaxies at 17 
Mpc (distance of the Virgo galaxies, dotted linos) , 
the resolution of 10 arcsec of Martin & Kennicutt 
(2001) correspond to 0.8 kpc, thus Vi is actually 
close to the quantity really involved in measure- 
ments. 

This star formation rate ipi at the threshold cu- 
riously shows very low scatter around 4 10^ Mq 
Gyr^^. The right axis of Fig. 8 shows V'l in the 
unit of number of massive stars present at any time 
within the 1 kpc wide annular ellipse. To compute 
this number, we integrated the initial mass func- 
tion (assumed to be a Salpeter function over the 
range 0.1 and 100 Mq) above 10 Mq, and assumed 
a lifetime for these stars (10^ yr). Our adopted 
limit of 10 Mq being conservative, our estimate 
should be an upper limit. Note that the predic- 
tion of the number of O stars using Starburst 99 
(Leitherer et al. 1999) with Solar mctallicity, and a 
Salpeter initial mass function between 1 and 100 
Mq is very close to the numbers we obtained in 
that way. Introducing lower mass stars into the 
initial mass function would reduce this number; 
thus once again, our number is probably an upper 
limit to the number of ionizing stars. 

Prom this Figure 8, we conclude that the 
threshold of star formation found by Martin & 
Kennicutt (2001) corresponds to the radius where 
the number of ionizing stars at any given time 
within 1 kpc wide ellipses becomes lower than 



about 1. At the position of the threshold, we 
would expect to find more than one star in only 
about 19 % of the galaxies of the galaxies of Martin 
& Kennicutt (2001). For three stars, the percent- 
age drops to about 6 %. It is thus natural that 
this radius will be the one where we stop finding 
HII regions, not because of the absence of star for- 
mation, but because there is a diminishingly small 
chance of catching a massive star in our annuli, 
even with a normal initial mass function. 

Although the threshold radius is defined by an 
abrupt drop in the Ha surface brightness pro- 
file, Martin & Kennicutt (2001) themselves no- 
ticed that many galaxies exhibit a few HII regions 
well beyond this radius, where star formation does 
seem more stochastic (see also the HII regions 
found in extreme outer parts of disk galaxies by 
Ferguson et al. 1998). The presence of resolved 
young blue B stars in the outskirts of M31 (Cuil- 
landre et al. 2001) also shows that some star for- 
mation occurs below the critical gas surface den- 
sity. Similarly, Davidge (2003) has shown that in- 
termediate age stars are presents at large radii in 
M33 and NGC2403, suggesting again that star for- 
mation has proceeded in the past beyond the cur- 
rent threshold radius. All these facts are totally 
consistent with our hypothesis that the threshold 
radius is only marking the last radius where HII 
regions are found in sufficient numbers to compute 
an Ha profile, and not the end of star formation. 

This is certainly relevant to the phenomenon 
of extended UV (XUV hereafter) disks in which 
UV emission is observed at large radii (and asso- 
ciated with relatively recent star formation) but 
where Ha detections are rare (Gil do Paz et al. 
2005; Thilker ot al. 2005). The situation is how- 
ever certainly not that simple since XUV disks are 
themsolvos not obvious extrapolations of the inner 
UV profiles. It will, of course, bo interesting to 
compare our results with the relation between gas 
and star formation rate in XUV galaxies. We plan 
to start that analysis with the spectacular case of 
NGC 4625 (Gil de Paz et al., in preparation). 

5. Conclusion 

Using the galaxies of large angular size culled 
from the GALEX Atlas of Nearby Galaxies (Gil 
de Paz et al. 2006), we have compiled GALEX 
UV and IRAS IR images to compute profiles of 
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attenuation and star formation in 43 galaxies. We 
find typically that UV attenuation decreases, from 
a few magnitudes of extinction in the centers of 
galaxies, to low values at larger radii. Combined 
with the metallicity gradients observed in those 
same galaxies, this gives rise to a attenuation- 
metallicity relationship, quantified in section 3.4. 
No clear correlation is found between the attenu- 
ation and the gas surface density. 

Using the attenuation to correct the FUV pro- 
files for dust, we derive star formation rate pro- 
files with which we test several modern expressions 
of the classical "Schmidt law" . Our data overlap 
with those of Kennicutt (1998a) but extend to sig- 
nificantly lower gas densities suggesting that the 
UV light in general traces and reveals lower levels 
of star formation than Ha, in coincidence with the 
recent GALEX discoveries of UV extended (XUV) 
disks (e.g. Gil de Paz et al. 2005; Thilker et al. 
2005). We suggest that the much-debated thresh- 
old of star formation might have been an observa- 
tional selection effect caused by the extremely low 
number of ionizing stars expected to be found be- 
yond the "threshold" radius. Indeed, for low levels 
of star formation, the Ha radiation in one ellipti- 
cal annuli seems stochastic when the number of 
massive stars in this annuli becomes close to ~ 1. 
The UV radiation coming from less massive and 
more long-lived stars (thus in larger number) does 
not present a break in its radial profile at the Ha 
threshold radius. 

While these results arc already interesting in 
their own right, it is obviously time to look to more 
galaxies, using Spitzer data in combination with 
GALEX. Spitzer will allow us to include many 
more galaxies in such studies owing to its bet- 
ter resolution than IRAS. Indeed, a joint analysis 
of GALEX and Spitzer data is already underway 
in a SINGS-GALEX collaboration. Note that the 
SINGS sample will only partially overlap with the 
one presented here as some of the closest galaxies 
were not included in SINGS. 

GALEX (Galaxy Evolution Explorer) is a 
NASA Small Explorer, launched in April 2003. 
We gratefully acknowledge NASA's support for 
construction, operation, and science analysis for 
the GALEX mission, developed in cooperation 
with the Centre National d'Etudes Spatiales of 
Prance and the Korean Ministry of Science and 



Technology. AGdP is partially financed by the 
Spanish Programa Nacional de Astronomfa y As- 
troffsica under grant AYA2003-01676. We also 
thank the MAGPOP network for its support, and 
the referee for very constructive comments. 
Facilities: GALEX 
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A. Appendix: individual profiles 



The loft panels of Figure 9 shows the profiles of UV attenuation, UV surface brightness (corrected for 
attenuation and inclination) equivalent to a star formation rate surface density, and of gas surface density 
used in this work. The right panels show gaseous hydrogen surface densities (HI, H2, HI+H2) vs the star 

formation rate surface density. 

[included in electronic version] 



18 



REFERENCES 

Aumann, H. H., Fowler, J. W., & Melnyk, M. 

1990, AJ, 99, 1674 

Bajaja, E., van der Burg, G., Faber, S. M., Gal- 
lagher, J. S., Knapp, G. R., & Shane, W. W. 
1984, A&A, 141, 309 

Bajaja, E., Wielebinski, R., Reuter, H.-P., Har- 
nett, J. I., & Hummel, E. 1995, A&AS, 114, 

147 

Barnes, J. E. 2004, MNRAS, 350, 798 

Bell, E. F. 2002, ApJ, 577, 150 

BeU, E. F., Gordon, K. D., Kcmiicutt, R. C., & 
Zaritsky, D. 2002, ApJ, 565, 994 

Belley, J., & Roy J.-R. 1992, ApJS, 78, 61 

Begeman, K. G. 1987, Ph.D. Thesis, 

Begeman, K. G. 1989, A&A, 223, 47 

Bersier, D., Blccha, A., Golay, M., & Martinet, L. 

1994, A&A, 286, 37 

Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, 
ApJ, 224, 132 

Boissier, S., & Prantzos, N. 1999, MNRAS, 307, 
857 

Boissier, S., & Prantzos, N. 2000, MNRAS, 312, 
398 

Boissier, S., Prantzos, N., Boselli, A., & Gavazzi, 
G. 2003, MNRAS, 346, 1215 

Boissier, S., Boselli, A., Buat, V., Donas, J., & 
Milliard, B. 2004, A&A, 424, 465 

Boissier, S., et al. 2005, ApJ, 619, L83 

Boselli, A., Gavazzi, G., Lequeux, J., Buat, V., 
Casoh, F., Dickey, J., & Donas, J. 1995, A&A, 
300, L13 

Boselli, A., Lequeux, J., & Gavazzi, G. 2002, 

A&A, 384, 33 

Boselli, A., Boissier, S., Cortese, L., Gil de Paz, A., 
Seibert, M., Madore, B.F., Buat, V., Martin, 
D.C., 2006, submitted to ApJ 

Bosma, A. 1978, Ph.D. Thesis, 



Bosma, A. 1981, AJ, 86, 1791 

Brocils, A. H., & van Woerden, H. 1994, A&AS, 
107, 129 

Braun, R., Walterbos, R. A. M., Kcnnicutt, R. C., 
& Tacconi, L. J. 1994, ApJ, 420, 558 

Braun, R. 1997, ApJ, 484, 637 

Buat, V. 1989a, A&A, 220, 49 

Buat, v., Dcharveng, J. M., & Donas, J. 1989b, 

A&A, 223, 42 

Buat, v., & Xu, C. 1996, A&A, 306, 61 

Buat, v., Boselh, A., Gavazzi, G., & Bonfanti, C. 
2002, A&A, 383, 801 

Buat et al., 2005 TO CHECK 

Burgarella, D., Buat, V., & Iglesias-Paramo, J. 
2005, MNRAS, 360, 1413 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, 
T. 1994, ApJ, 429, 582 

Calzetti, D., et al. 2005, ApJ, 633, 871 

Carignan, C, & Puche, D. 1990, AJ, 100, 394 

Carignan, C, & Puche, D. 1990, AJ, 100, 641 

Casertano, S., & van Gorkom, J. H. 1991, AJ, 101, 
1231 

Cayatte, V., Kotanyi, C, Balkowski, C, & van 
Gorkom, J. H. 1994, AJ, 107, 1003 

CorbeUi, E. 2003, MNRAS, 342, 199 

Cortese, L., et al. 2006, ApJ, accepted, astro- 
ph/0510165 

Cuillandre, J.-C, Lequeux, J., Allen, R. J., Mel- 
lier, Y., & Bertin, E. 2001, ApJ, 554, 190 

Dale, D. A., Helou, G., Contursi, A., Silbermann, 
N. A., & Kolhatkar, S. 2001, ApJ, 549, 215 

Dame, T. M., Koper, E., Israel, F. P., & Thad- 
deus, P. 1993, ApJ, 418, 730 

Davidge, T. J. 2003, AJ, 125, 3046 

Donas, J., Deharveng, J. M., Laget, M., Milliard, 
B., & Huguenin, D. 1987, A&A, 180, 12 



19 



Dopita, M. A., & Ryder, S. D. 1994, ApJ, 430, 163 

Elmegreen, B. G. 2002, ApJ, 577, 206 

Elmegreen, G., & Hunter, D., 2006, to appear in 
ApJ, astro-ph/0509190 

Ferguson, A. M. N., Wyse, R. F. G., Gallagher, 
J. S., & Hunter, D. A. 1998, ApJ, 506, L19 

Gavazzi, G., Zaccardo, A., Sanvito, G., Boselli, 
A., & Bonfanti, C. 2004, A&A, 417, 499 

Gil de Paz, A., et al. 2005, ApJ, 627, L29 

Gil dc Paz ct al., 2006, in press, astro-pli/0606440 

Guelin, M., Zylka, R., Mezger, P. G., Haslam, 
C. G. T., Kreysa, E., Lemke, R., & Sievers, 
A. W. 1993, A&A, 279, L37 

Hatton, S., Devriendt, J. E. G., Ninin, S., 
Bouchct, F. R., Guiderdoni, B., & Vibert, D. 

2003, MNRAS, 343, 75 

Heckman, T., et al. 1995, ApJ, 452, 549 

Heckman, T. M., Robert, C, Leitherer, C, Gar- 
nett, D. R., & van der Rydt, F. 1998, ApJ, 503, 
646 

Heyer, M. H., Corbelli, E., Schneider, S. E., & 
Young, J. S. 2004, ApJ, 602, 723 

Holwcrda, B. W., Gonzalez, R. A., Allen, R. J., & 
van der Kruit, P. C. 2005, AJ, 129, 1381 

Holwerda, B. W., Gonzalez, R. A., Allen, R. J., & 
van der Kruit, P. C. 2005, A&A, 444, 101 

Hunter, D. A., Elmegreen, B. G., & van Woerden, 
H. 2001, ApJ, 556, 773 

Jackson, D. C., Skillman, E. D., Cannon, J. M., & 
Cote, S. 2004, AJ, 128, 1219 

Jansen, R. A., Fabricant, D., Franx, M., & Cald- 
weh, N. 2000, ApJS, 126, 331 

Kennicutt, R. C. 1989, ApJ, 344, 685 

Kennicutt, R. C. 1998a, ApJ, 498, 541 

Kennicutt, R. C. 1998b, ARA&A, 36, 189 

Kennicutt, R. C, Bresohn, F., & Garnett, D. R. 
2003a, ApJ, 591, 801 



Kennicutt, R. C, et al. 2003b, PASP, 115, 928 

Kenney, J. D., & Young, J. S. 1988, ApJS, 66, 261 

Kenney, J. D. P., Scoville, N. Z., & Wilson, C. D. 
1991, ApJ, 366, 432 

Komugi, S., Sofuc, Y., Nakanishi, H., Onodera, S., 
& Egusa, F. 2005, PASJ, 57, 733 

Kong, X., Chariot, S., Brinchmann, J., & Fall, 
S. M. 2004, MNRAS, 349, 769 

Lake, G., & Skillman, E. D. 1989, AJ, 98, 1274 

Lee, H., Grebel, E. K., & Hodge, P. W. 2003, 
A&A, 401, 141 

Lee, H., McCall, M. L., Kingsburgh, R. L., Ross, 
R., & Stevenson, C. C. 2003, AJ, 125, 146 

Leitherer, C, et al. 1999, ApJS, 123, 3 

Lia, C, Portinari, L., & Carraro, G. 2002, MN- 
RAS, 330, 821 

Madore, B. F. 1977, MNRAS, 178, 1 

Martimbeau, N., Carignan, C, & Roy, J.-R. 1994, 
AJ, 107, 543 

Martin, C. L., & Kennicutt, R. C. 2001, ApJ, 555, 
301 

Masegosa, J., Moles, M., & del Olmo, A. 1991, 
A&A, 249, 505 

Meurer, G. R., Heckman, T. M., Leitherer, C, 
Kinney, A., Robert, C, & Garnett, D. R. 1995, 

AJ, 110, 2665 

Meurer, G. R., Heckman, T. M., & Calzetti, D. 
1999, ApJ, 521, 64 

Morrissey, P., et al. 2005, ApJ, 619, L7 

Mulder, P. S., & van Driel, W. 1993, A&A, 272, 
63 

Newton, K., & Emerson, D. T. 1977, MNRAS, 
181, 573 

Nictcn, C, Neininger, N., Guelin, M., Ungerechts, 
H., Lucas, R., Berkhuijsen E.m., Beck R., 
Wielebinski R. 2006, A&A, 453, 459 

Ohnishi, T. 1975, Progress of Theoretical Physics, 
53, 1042 



20 



Ondrcchcn, M. P., van dcr Hulst, J. M., & Hum- 
mel, E. 1989, ApJ, 342, 39 

Paglione, T. A. D., et al. 2001, ApJS, 135, 183 

Panuzzo, P., Bressan, A., Granato, G. L., Silva, 
L., & Dancsc, L. 2003, A&A, 409, 99 

Pisano, D. J., Wilcots, E. M., & Elmegreen, B. G. 
1998, AJ, 115, 975 

Popescu, C. C, Misiriotis, A., Kylafis, N. D., 
Tuffs, R. J., & Fischera, J. 2000, A&A, 362, 

138 

Prantzos, N., & Aubert, O. 1995, A&A, 302, 69 

Puche, D., Carignan, C, & Bosma, A. 1990, AJ, 
100, 1468 

Puche, D., Carignan, C., & van Gorkom, J. H. 
1991, AJ, 101, 456 

Puche, D., Carignan, C, & Wainscoat, R. J. 1991, 
AJ, 101, 447 

Quirk, W. J. 1972, ApJ, 176, L9 

Rand, R. J. 1994, A&A, 285, 833 

Regan, M. W., Thornley, M. D., Heifer, T. T., 
Sheth, K., Wong, T., Vogcl, S. N., Blitz, L., & 
Bock, D. C.-J. 2001, ApJ, 561, 218 

Rice, W. 1993, AJ, 105, 67 

Roy, J.-R., Belley, J., Dutil, Y., & Martin, P. 1996, 
ApJ, 460, 284 

Sakamoto, S., Handa, T., Sofue, Y., Honma, M., 
& Sorai, K. 1997, ApJ, 475, 134 

Sancisi, R., & Allen, R. J. 1979, A&A, 74, 73 

Schaye, J. 2004, ApJ, 609, 667 

Schiminovich, D., et al. 2005, ApJ, 619, L47 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 
1998, ApJ, 500, 525 

Schmidt, M. 1959, ApJ, 129, 243 

Seibert, M., et al. 2005, ApJ, 619, L55 

Seigar, M. S. 2005, MNRAS, 361, L20 

Shostak, G. S., & van der Kruit, P. C. 1984, A&A, 
132, 20 



Storchi-Bergmann, T., Wilson, A. S., & Baldwin, 
J. A. 1996, ApJ, 460, 252 

Sofue, Y. 1997, PASJ, 49, 17 

Sorai, K., Nakai, N., Kuno, N., Nishiyama, K., & 
Hascgawa, T. 2000, PASJ, 52, 785 

Thilker, D. A., et al. 2005, ApJ, 619, L79 

Thilker et al., 2006, in preparation 

Toomre, A. 1964, ApJ, 139, 1217 

van Dricl, W., Rots, A. H., & van Woerden, H. 

1988, A&A, 204, 39 

van Driel, W., Valentijn, E. A., Wessclius, P. R., 
& Kussendrager, D. 1995, A&A, 298, L4 Van 
Driel et al. 95 FOR N0660 ; 

van Zee, L., Salzer, J. J., Haynes, M. P., 
O'Donoghue, A. A., & Balonek, T. J. 1998, AJ, 

116, 2805 

Vollmer, B., & Bcckcrt, T. 2002, A&A, 382, 872 
Wang, B., & Silk, J. 1994, ApJ, 427, 759 

Wevcrs, B. M. H. R., van dcr Kruit, P. C, & Allen, 

R. J. 1986, A&AS, 66, 505 

Wilding, T., Alexander, P., & Green, D. A. 1993, 
MNRAS, 263, 1075 

Witt, A. N., & Gordon, K. D. 2000, ApJ, 528, 799 

Wong, T., & Blitz, L. 2002, ApJ, 569, 157 

Wyse, R. F. G., & Silk, J. 1989, ApJ, 339, 700 

Young, J. S., & ScoviUe, N. 1982, ApJ, 260, L41 

Young, J. S., Tacconi, L. J., & Scoville, N. Z. 1983, 
ApJ, 269, 136 

Young, J. S., et al. 1995, ApJS, 98, 219 

Zaritsky, D., Kennicutt, R. C, & Huchra, J. P. 
1994, ApJ, 420, 87 

Zhang, X., Wright, M., & Alexander, P. 1993, 
ApJ, 418, 100 



This 2-column preprint was prepared witii the AAS lAT^^X 
macros v5.2. 



21 



-3-2-1012 -3-2-1012 




-0.5 0.5 1 1.5 2 -0.5 0.5 1 1.5 2 

FUV-NUV FUV-NUV 

Fig. 1. — Total infrared to FUV ratio vs FUV-NUV color. The FUV-NUV color is a measure of the slope in 
the UV f3 given at the top (Kong et al. 2004). The right axis indicates the scale in A{FUV) using equation 
3. In the left panel, points with errorbars larger than 0.5 (dex on left axis and mag on bottom axis) were 
omitted for the sake of clarity (errorbars indicate how much the results would change if our sky determination 
in the UV or far infrared was moved by ±lcr). The line is the usual relation for startbusts as given in Kong et 
al. (2004). (Electronic version: different colors indicates various axis ratio b/a: red for < b/a < 1/3; black 
for 1/3 < b/a < 2/3; blue for 2/3 < b/a < 1). Grey squares are the integrated values of the GALEX Atlas 
of nearby galaxies (Gil de Paz et al. 2006). Squares with a circle correspond to the galaxies in our study. 
In the right panel, we show as black squares the points used to compute our fits (excluding nearly edge-on 
galaxies). The usual starbursts relationship is compared to our fit and its 1 sigma deviation for a simple 
fit (dashed curve and horizontally hatched area) and a fit taking into account the uncertainties (dotted line 
and vertically hatched area). 
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Fig. 2. — Attenuation at FUV wavelength {FtibJ Ffuv ratio scale on the right axis) as a function of the 
observed column density of HI (left), H2 (right). Points belonging to a same galaxy are connected by a 
dotted line. Different symbols indicate different values of the axis ratio (b/a) as indicated in the corner. 
Note that since A{FUV) is not corrected for inclination, the gas column densities used correspond to the 
same inclination (hence the "inc" subscript). Errorbars indicate how much the results would change if our 
sky determination in the UV or far infrared was moved by ±1(T. 
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Fig. 3. — Attenuation at FUV wavelength 
{FtirJ Fpuv ratio scale on the right axis) as 
a function of the observed total gas density 
(HI+H2). Since A{FUV) is not corrected for 
inclination, the gas column densities used corre- 
spond to the same inclination (hence the "inc" 
subscript). Errorbars indicate how much the re- 
sults would change if our sky determination in the 
UV or far infrared was moved by ztltr. 
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Fig. 4. — Attenuation at FUV wavelength 
{Ftir/ Fpuv ratio scale on the right axis) as a 
function of the metallicity (as given by the abun- 
dance gradient of each galaxy). The short-dashed 
line is a fit (A(FUV)=1.02 (12-Hlog(0/H)) - 7.84) 
to the points belonging to galaxies with b/a > 0.4 
(see section 3.4). The long-dashed curve corre- 
sponds to the fit ol FtirJ Fpuv vs 0/H in Cortese 
et al. (2006). Errorbars indicate how much the re- 
sults would change if our sky determination in the 
UV or far infrared was moved by ±lcr. 
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Fig. 5. — Test of various "Schmidt laws" for our 
whole sample. The star formation rate (right axis) 
derived from FUV profiles corrected for attenua- 
tion and inclination is compared to HI (top left), 
H2 (top right), total H surface densities (bottom- 
left). In the bottom- right panel, we show the SFR 
X R/V(R) vs hydrogen total density to test equa- 
tion 6. In each box, the hourglass shaped shaded 
area indicates equivalent results from Boissier et 
al. (2003) based on Ha profiles (the extent corre- 
sponds roughly to the extent of the observations 
and the diagonals to the regression lines of this 
study). The correlation coefficient and regression 
lines (dashed lines, obtained fitting y vs x and x 
vs y, see section 4.2) are shown (their parameters 
are given in table 2). Errorbars indicate how much 
the results would change if our sky determination 
in the UV or far infrared was moved by ilcr. 



Fig. 6. — The simple Schmidt law (Star forma- 
tion rate surface density vs total Hydrogen sur- 
face density) found in our work along radial pro- 
files of GALEX galaxies (including the extinction- 
correction: open circles, or not: crosses) is com- 
pared to the one of Kennicutt (1998a) (filled cir- 
cles, obtained by computing average surface den- 
sities within the optical disk). The shaded area 
indicates gas densities values found at the thresh- 
old by Kennicutt (1989). The studies overlap with 
each other in the high density half, but we obtain 
many more points at low surface densities. The 
points corresponding to M31 (following an untyp- 
ical path in the inner galaxy) are connected: high 
SFR and low gas amounts (top-left) correspond to 
the inner galaxy. When increasing the radius, the 
SFR decreases while the gas increases (contrary 
to classical SFR "laws"). At about 50 arcmin, the 
trend changes and both the SFR and the gas sur- 
face densities go down with radius, following the 
SFR "law" described by other galaxies. 
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Fig. 7. — Profiles of star formation rate surface density determined from the UV and corrected for extinction 
with the far infrared (dotted) and from Ha (solid, taken from Martin and Kennicutt, 2001) for the 9 galaxies 
in common with Martin and Kennicutt (2001), having a threshold radius larger than 90 arcsec (so that 
our data resolves it). The radius is normalised by the "threshold radius" Rthreshoid measured by Martin 
and Kennicutt (2001). We show only a region around the threshold radius, marked by a vertical dashed 
line. Note that for the Ha profiles, the horizontal parts of the profiles beyond R/ Rthreshoid ^ 1 are actually 
upper limits. The Ha profiles have a much higher resolution (10 arcsecs) than the UV + far infrared data 
(resolution shown by the errorbar in the central-bottom part of each panel), however we have many UV 
points beyond the threshold without any sign of a decrease in the UV for most galaxies. 
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Fig. 8. — Profiles of tlie star formation rate surface 
density determined from the Ha profiles of Martin 
& Kennicutt (2001) (MKOl). The left axis indi- 
cates the star formation rate integrated in 1 kpc 
annular ellipses. The right axis shows the num- 
ber ionizing stars within the same annular ellipses 
(see text for details). Dotted lines correspond to 
galaxies in Virgo. The radius is normalised by the 
"threshold radius" Rthreshoid measured by Mar- 
tin & Kennicutt (2001). We show only a region 
around the threshold radius, marked by a verti- 
cal dashed line. Note that several points beyond 
R/ Rthreshoid (straight almost horizontal lines) are 
actually upper limits. 



27 



> 

El, 



22 

^ 24 

(Q 

J, 

5 28 

D 

* 30 
32 

8 

a) 

3 

: 6 
4 
2 


10 



© 



0.1 



0.01 



; ■■■■■-■■} 


i 1 1 1 1 i 1 1 
1 1 ^ 1 1 1 1 1 i 


' 1 1 1 1 1 1 1 1 1 1 1 1 1 

Xi si .Zi iSl 




total 

HI 

' . 

... 1 ... 1 ... 1 ... 1 ... 1 . . 



200 



400 600 
R (") 



800 1000 




0.1 1 



Fig. 9.— [THIS IS ONE PANEL OF FIG. 
9. THE COMPLETE FIG. 9 WILL BE 
AVAILABLE IN ELECTRONIC FORM 
ON THE APJ WEB SITE AND AT: 
http: / /www. oamp.fr/boissier/preprint / prcprint.html 
] Radial profiles of the gas surface density (bottom- 
left), FUV attenuation (middle-left), FUV surface 
brightness (corrected for attenuation and to a 
face-on inclination, top-left) when available. This 
FUV surface brightness is directly proportional 
to the star formation rate surface density shown 
in the right panels as a function of the HI, H2 
and the total gasesous Hydrogen surface density 
along those profiles. Errorbars in the left panels 
indicate how much the results would change if 
our sky determination in the UV or far infrared 
was moved by iltr (see section 2). The name of 
the galaxy is indicated at the top. 
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